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We present lattice QCD simulation results from the European Twisted Mass Collabo- 
ration (ETMC) for the light, strange and charm quark contents of the nucleon. These 
quantities are important ingredients to estimate the cross-section for the detection of 
WIMPs as Dark Matter candidates. By employing a particular lattice QCD formulation, 
i.e. twisted mass fermions, accurate results of the light and strange scalar contents of the 
nucleon can be obtained. In addition, we provide a bound for the charm quark content 
of the nucleon. 
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1. Introduction 

Dark Matter is still quite mysterious and its nature not clarified. Quite popular 
candidates are weakly interacting massive particles (WIMPs) which appear in many 
models pQ. There are a number of experiments that attempt to detect such WIMPs. 
The underlying mechanism for the WIMP detection is that -due to their assumed 
large mass- they produce a Higgs boson which in turn can interact with nucleons. 
What the experiments are -in principle- able to detect then is the recoil energy of 
the nucleon struck in this way by a WIMP which is flying through the detector. 

Since the Higgs boson is a scalar, the proposed interaction of the WIMP and the 
nucleon is going to occur via the coupling of the Higgs boson to the scalar quark 
content of the nucleon. The spin independent (SI) elastic cross section for such a 
process reads [5] 



where G q . is the effective coupling constant between a quark of flavour / and 
the WIMP of mass m x and {N\qjqf\N) denotes the scalar content of the nucleon 
belonging to a quark flavour /. The strength of the coupling G q . will depend on 
the particular model it is computed in. 

It is important to note that changes of the scalar content of the nucleon by 
about 10% can lead to large changes in the cross-section since (N\q~fqf\N) 

appears quadratically in eq. ([T]). One problem with calculating the cross-section is 
that the scalar quark contents of the nucleon are only poorly determined. If chiral 
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perturbation theory (%PT) is used for their calculations one hnds for the ratio tjn 
of the strange to the light scalar contents 

2{N\ss\N) 

VN = (N\uu + dd\N) ( ' 

a value of y N = 0.44(13) 0]. This suggests that actually the strange quark content 
of the nucleon is surprisingly large. On the other hand, from this X-PT evaluation of 
Un only an accuracy of about 30% is achieved leading to a large uncertainty in the 
cross-section itself. Thus, is seems highly desirable to compute yN accurately and 
test whether indeed the scalar quark content of the nucleon is as large as indicated 
by xPT. 

In principle, lattice QCD provides a tool to reach a precise quantitative result 
for ypj. However, also in this framework the computation of yN turns out to be very 
difficult. One way to calculate yN is via the Feynman-Hellman theorem |o which 
uses the quark mass dependence of the nucleon mass, 

ovjv = m(N\uu + dd\N) = m N and <j ss n — m s (N\ss\N) — m s — — — , (3) 

kJ I I L (J II ig 

where we denote with m the mass-degenerate light up and down quark masses and 
with m s the strange quark mass. However, to determine yN taking the numerical 
derivative requires many simulations at different values of the quark masses. In 
addition, it is unclear how sensitive the nucleon mass will be on the variation of the 
quark mass. 

A second approach is to compute the scalar quark contents (N\q~fqf\N) directly. 
The disadvantage here is that dis-connected (singlet) contributions are to be com- 
puted which are generically very noisy in lattice calculations and need a very large 
statistics. 

In this proceedings contribution we want to report about a lattice determination 
of yN which uses a particular formulation of lattice QCD, namely the so-called 
twisted mass fermions taken at maximal twist. What is important for the discussion 
here is first of all that physical quantities computed within this approach scale to the 
continuum limit, i.e. sending the lattice spacing a to zero, with a rate that is of 0(a 2 ) 
leading to a rapid approach towards the continuum limit. Secondly, twisted mass 
fermions allow for special techniques to compute dis-connected diagrams which give 
a much better signal to noise ratio than most other lattice formulations of QCD. 
In this report we will make heavily use of this special property of twisted mass 
fermions. A third advantage of twisted mass fermions (at maximal twist) is that 
they allow for a straightforward renormalization of (N\qfqf\N) [6] avoiding any 
mixing between the bare light, strange and charm quark matrix elements. Although 
also chiral invariant lattice fermions share the same property, these latter kind of 
lattice fermions are much more computer time demanding such that twisted mass 
fermions, from this point of view, are clearly advantageous. 

We want to mention that there exist already a number of computations of the 
strange quark content of the nucleon, see refs. [7l-ll4j. For a review we refer to ref. [4]. 
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However, because of the aforementioned reasons, these calculations are affected with 
rather large statistical errors. It is our goal to improve significantly on this situation 
and provide a more precise value of the scalar content of the nucleon for the light 
and strange quarks. In addition, we will give a bound on the charm quark content. 

2. The calculation and the results 

As said above, we use a particular formulation of lattice QCD called twisted 
mass fermions as proposed by Frezzotti and Rossi, see refs. [T5HT5] . The Euro- 
pean Twisted Mass Collaboration (ETMC) -in which this project is embedded- is 
using these kind of lattice fermions for a number of years now and has obtained 
many physical results already, see e.g. refs. [T9H24] and references therein. 

In the latest set of simulations, ETMC has incorporated the complete first two 
quark generations in the simulations, meaning that besides the light up and down 
quarks also the strange and charm quarks are taken into account as dynamical 
degrees of freedom |24H26| . This is presently worldwide a rather unique setup for 
lattice simulations and particularly important for the problem addressed here since 
we want to evaluate the strange and the charm quark contents of the nucleon. 
In order to compute nucleon scalar matrix elements a nucleon 2-point function is 
needed which is given by 

Cw,2pts(* - tsrc) = tJ± ( Jn(x)Jn~(x stc )), (4) 

x 

where the subscript N refers to the proton or to the neutron states for which the 
interpolating fields are given by : 

JV = e abc ( u a,T Cl5d b^ ^ ^ jn = £ a6c (^,^75^) d C . 

The projectors used are T ± = 1= ^ 7 ° , and C is the charge conjugation matrix. Using 
discrete symmetries and anti-periodic boundary conditions in the time direction for 
the quark fields, we have 2pts (i) = —Cjj 2 pts(^ — 

The other input for the computation of the scalar quark content of the nucleon 
is the zero momentum 3-point function 

Cw;St B (*«» A *qp) - E ^ ± (JN(x)O q (x op )J N '(x src )), (5) 

where O q is an operator having scalar quantum numbers, At op = t op — t SIC is 
the time of insertion of the operator, and t s — t — t src gives the so-called source- 
sink separation. Since we will consider an operator with a non vanishing vacuum 
expectation value, we also define 

CjV,'3pts A *°p) = CjV,'3pts(^' At op ) - C^,2pts(*' X src) ^(Qg^op)) • (6) 

The scalar quark operators O q are given by 

Oi = uu + dd, O s — ss or O c = cc, (7) 
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depending on the quantity of interest. Using the definitions of the two- and three- 
point functions of Eqs. ((4]) and ([5]), we will consider the ratio 

r ±,O q ,vcv, s 
K0 q {ts,t O p) - — q: — (») 

which is related to i/n and will serve as our test quantity for the noise reduction 
technique to compute dis-connected contributions. Our analysis has been performed 
on a 32 3 x 64 volume with a lattice spacing of a = 0.0777(4) fm and a pion mass 
of approximately 380 MeV. In order to improve the overlap between the ground 
state and the interpolating operators we use Gaussian smearing of the quark fields 
appearing in the interpolating fields. We use APE smearing of the gauge links 
involved in the Gaussian smearing, following the same strategy as in |2"71l2"8] . 

We will use the powerful variance reduction method for twisted mass fermions 
introduced in [29,30 and used to study the 7/ meson in [31]. In Fig.[T] we compare 
the efficiency of this technique with another, more standard noise reduction tech- 
nique which relies on the hopping parameter expansion of the Dirac operator. This 
latter technique can be used for any discretization, and has been introduced in [32] ■ 
We refer the interested reader to appendix B. of ref. [30] for the implementation of 
the method for twisted mass fermions. The here used new variance noise reduction 
technique improves the signal to noise ratio by a factor ~ 3 after having fixed all 
other parameters. In particular, as can be seen in fig. [T]with a statistics of 0(1000) 
gauge field configurations we can reach a result that is at the 5tr level different from 
zero, while the standard method with such a statistics gives only a result at the la 
level. 

Using our improved noise reduction technique for dis-connected contributions, 
we have first performed an investigation of the light quark content of the nucleon, 
i.e. ovat of eq. ([3]). Here, it turns out that the connected part of the relevant 3- 
point function are by far dominating and that the dis-connected contribution can 
be safely neglected. Although such a fact might have been expected, our dedicated 
investigation is certainly important to confirm this. 

To estimate systematic effects coming from the excited states, the final number 
for ovat is then obtained from a distribution of many fit results. From the (weighted) 
mean and variance we then obtain 

ovaK^ps ~ 380MeV) = 150(l)(10)MeV , (9) 

where the first error is statistical and the second gives the systematic uncertainties. 
We then proceeded to compute the ratio R of eq. (|8|) for the strange and charm 
quark scalar operator. In Fig. [2] (left), we show the plateau corresponding to the 
strange quark. The plateau value appears to be 5.7 a away from zero showing that 
we indeed can obtain a rather accurate value for the strange quark content of the 
nucleon. 

We also carried out a first study of the charm quark content of the nucleon. Here 
it is most important that we can use gauge field configurations from Nj = 2 + 1 + 1 
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Fig. 1. Signal to noise ratio of the correlator ratio Eq. (§). as a function of the number of 
configurations used, TV, for the variance noise reduction technique used in this paper compared 
to the more standard hopping parameter expansion technique. Note that here the source-sink 
separation is fixed to t a = 12a. Curves are only shown to guide the eyes. Our method allows to 
reach results at a ~ 5<r significance level with a statistics that is O(10) less than standard methods. 
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Fig. 2. In the left graph, we plot the bare ratio of Eq. {8} for the strange quark. We used 842 
configurations for this analysis. Note that the signal is ha away from 0. In the right graph, we 
show results for the charm content using the same statistics. While certainly no plateau can be 
identified, the charm quark content is clearly not larger than the strange quark content of the 
nucleon. As in fig. [T] t s = 12a. 



simulations that take into account also the charm quark degree of freedom. We show 
in Fig. [2] (right) our results for the charm quark. Unfortunately, there is no plateau 
signal visible. The gray band in Fig. [5] (right) represents our results for the strange 
quark. The graph then suggests that the charm quark content of the nucleon cannot 
be large compared to the strange quark as one might expect. 
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3. Conclusion 

In this proceedings contribution we have computed the scalar quark content of the 
nucleon for all four quarks of the first two generations. We find a value of the 
phenomenologically important quantity i/n = 0.066(H)(2), see eq. ([2]). This value 
is significantly smaller than predicted in chiral perturbation theory and indicates 
that the strange quark content is only a few percent. In fig.[3]we give a summary of 
presently available values for y N demonstrating the precision of our result. We also 
could provide for the first time a bound for the charm quark content of the nucleon 
excluding that it is significantly larger than the strange quark content. 
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Fig. 3. Compilation of presently available results for the y-parameter jyjv of eq. J2}. 

As a present shortcoming, we mention that our calculations are performed at 
only one value of the lattice spacing, quark mass and volume. Thus we are lack- 
ing presently a good estimate of systematic effects originating from varying these 
parameters of the simulations. Clearly, we will address these systematic errors in 
future computations which have already been started. 
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